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Abstract 

Objective:lo project the HIV/AIDS epidemics among men who have sex with men (MSM) under different combinations of 
HIV testing and linkage to care (TLC) interventions including antiretrovlral therapy (ART) In Beijing, China. 

Design: Mathematical modeling. 

Methods: Us\r\g a mathematical model to fit prevalence estimates from 2000-2010, we projected trends in HIV prevalence 
and incidence during 2011-2020 under five scenarios: (SI) current intervention levels by averaging 2000-2010 coverage; 
(S2) increased ART coverage with current TLC; (S3) increased TLC/ART coverage; (S4) increased condom use; and (S5) 
increased TLC/ART plus increased condom use. 

Results: The basic reproduction number based upon the current level of interventions is significantly higher than 1 
(i^o =2.09; 95% confidence interval (CI), 1.83-2.35), suggesting that the HIV epidemic will continue to increase to 2020. 
Compared to the 2010 prevalence of 7.8%, the projected HIV prevalence in 2020 for the five prevention scenarios will be: 
(SI) Current coverage: 21.4% (95% CI, 9.9-31.7%); (S2) Increased ART: 19.9% (95% CI, 9.9-28.4%); (S3) Increased TLC/ART: 
14.5% (95% CI, 7.0-23.8%); (S4) Increased condom use: 13.0% (95% CI, 9.8-28.4%); and (S5) Increased TLC/ART and condom 
use: 8.7% (95% CI, 5.4-11.5%). HIV epidemic will continue to rise (_Ro>l) for S1-S4 even with hyperbolic coverage in the 
sensitivity analysis, and is expected to decline (i?o=0.93) for S5. 

Conclusion: Our transmission model suggests that Beijing MSM will have a rapidly rising HIV epidemic. Even enhanced 
levels of TLC/ART will not interrupt epidemic expansion, despite optimistic assumptions for coverage. Promoting condom 
use is a crucial component of combination interventions. 
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Introduction 

More than 30 years have passed since the first HIV case was 
reported among men who have sex with men (MSM) [1]. MSM 
contribute to the majority of new cases in tlie Americas, Western 
Europe, Oceania, and much of Asia [2]. Of tlie 780,000 people 
living with HIV/ AIDS in China, 17.4% were infected through 
homosexual contacts in 201 1, rising from 1 1.0% in 2007 and 7.3% 
in 2005 [3]. A recent survey among 47,231 MSM from 61 Chinese 



cities showed an overall prevalence of 4.9% [4]. The unrelenting 
HIV burden among MSM suggests that HIV intervention remains 
an urgent priority. One promising concept is the expansion of 
HIV testing and linkage to care (TLC) including antiretrovlral 
therapy (ART), so that a much higher proportion of infected men 
know their HIV status and undergo ART, thereby reducing their 
infectiousness to others [5,6,7,8,9]. 
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China has implemented several aggressive public health 
programs in the past decade, including routine HIV testing and 
risk reduction interventions for high risk populations and free 
ART for HIV-infected treatment-eligible individuals [10]. These 
programs curtailed HIV transmission through unhygienic plasma 
collection and injecting drug use, and reduced mortality among 
HIV-infected patients [11,12,13]; however, no programs have 
reversed the rising HIV epidemic among Chinese MSM [14,15]. 
A recently published mathematic model estimated that a four-fold 
increase in testing rates may prevent 42,000 new HIV infections 
among aU at-risk groups in China over 5 years [16]. A few 
mathematical models evaluated independent or joint impacts of 
HIV testing, risk reduction intervention and ART among MSM 
[17,18,19,20,21,22,23,24], but these studies were mainly from the 
Americas, Western Europe, and Australia, where the HIV 
epidemics have been long established among MSM, and none 
from regions with a rapidly rising HIV epidemic among MSM. 
ART may not have the same protection against HIV transmission 
in homosexual contacts as in heterosexual contacts due to 
biological and behavioral differences in infectiousness/ susceptibil- 
ity or differences in co-factor frequency [25]. Furthermore, the 
eflFectiveness of combined TLC intervention packages may be 
contingent on the epidemiological context [26]. 

In this study, we modeled the dynamics of the HIV/ AIDS 
epidemic among MSM in China's capital city, Beijing, which is 
representative of most Chinese urban areas [4]. We simulated 
HIV prevalence from 2000-2010 based on the coverage and 
intensity of existing public health programs including uptake of 
HIV testing, linkage to HIV care (e.g. risk reduction), ART and 
condom use. Further, we made projections for the HIV epidemic 
beyond 2010 under various intervention scenarios. We then 
examined the robustness of model predictions based on different 
underlying assumptions in sensitivity analyses. 

Methods 

We used a compartmental ordinary differential equations model 
for simulating and projecting the HIV epidemic among MSM in 
Beijing. Instead of the bilinear incidence [27,28,29] and the 
standard incidence models [5,8,9], we used Bernoulli processes to 
describe the probability of HIV transmission [30]. Bernoulli processes 
can more accurately describe the HIV transmission network 
among MSM by giving consideration for multiple sexual partners, 
type of sexual partner (i.e., regular or casual), differential condom 
use, and probability of transmission by preferred anal sex position 
(i.e., receptive or insertive). Substance abuse is rare among MSM 
in Beijing; therefore, injecting drug use is unlikely to contribute 
substantially to HIV transmission [31,32,33]. Additionally, some 
Chinese MSM may have female sexual partners, but these female 
partners generally do not have risky behavior and are unlikely to 
transmit disease to male partners [34]. Therefore, we only 
considered HIV transmission as a result of homosexual contact; 
we implicitly assumed there is no other route of transmission 
[31,33,35]. We also assumed that the age of active homosexual 
intercourse ranges from 18-60 years among Chinese MSM [36]. 

Model structure 

Receipt of HIV testing and risk reduction services may change 
the risk of HIV transmission among MSM. Thus, we divided the 
study population into three mutually exclusive sub-populations: 
those who are not testing for HIV (Non-testing), those who are 
testing but not reducing their risk of infection/transmission 
(Testing), and those who are testing and reducing their risk of 
infection/transmission (Risk reduction). The population was 



further compartmentalized by infection status and disease status 
(or eligibility for ART treatment). According to Chinese national 
ART guidehnes [37], patients with CD4-H T cell count <350 cells/ 
|a.L are eligible for ART treatment. 

Figure 1 shows the schematic diagram of our compartmental 
model structure. In this model we assumed that uninfected MSM 
enter the model as susceptible [Sm) (these MSM could be local 
Beijing residents who turn 18 years old, or those who migrate from 
other areas to Beijing). Some who receive HIV testing may test 
negative {St), and then a proportion of those HlV-negative men 
are linked to care and receive risk reduction intervention provided 
by public health programs in Beijing {Sr). The subscripts N,T,R 
refer to the non-testing, testing, and risk reduction, respectively. 

Upon infection, individuals move from susceptible into infected 
status, and the latter is further divided into two infection stages: 
CD4 count >350/|aL (4,4 and 4) and CD4 count <350/|.lL 
(4>4 ^^'^ ^r)- Flow is unidirectional from susceptible status to 
infection stage 1 and then 2, as without ART, an HIV-infected 
individual is unlikely to experience immunological recovery. A 
proportion of infected MSM (4 or 4) receive ART (4 or 4); 
these compartments are bidirectional due to attrition or adherence 
to ART and decline or improvement of immunological status. 

The key to quantifying HIV incidence among MSM in Beijing 
is specifying the time-dependent HIV transmission rate, /.(?). 
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Figure 1. Schematic diagram of HIV combination prevention 
intervention model in the presence of ART. The top row of 

compartments show uninfected IVISM who are Susceptible to HIV 
infection. They are divided into Non-testers, Testers, and testers who 
received Risk reduction and linkage-to-care interventions. The bottom 
two rows representing two stages of disease contain an additional 
compartment — Infected individuals. Stage 1 is CD44- cell count <350/ 
|j.L while stage 2 is CD44- >350/|j.L The last column of compartments 
shows HIV-infected individuals who have initiated Antiretroviral 
therapy, a is the recruitment rate per year into the MSIVl population 
in Beijing. They are local HIV- MSM who turn 18 years old, or those from 
other areas who immigrate to Beijing. Horizontal movement is 
parameterized by rates of testing (k), rates of linkage to care (/), rates 
of ART initiation (i), and rates of ART failure or dropout (i//). Movement 
from Susceptible to Infected U) is the HIV transmission rate. It is 
estimated from multiple parameters including condom use, insertive/ 
receptive intercourse, regular/casual partner encounter, CD4 count, and 
reductions in infectivity due to TLC or ART. The natural removal rate (^i) 
represents death, being older than 60 years, or migration out of Beijing. 
The rate of disease progression is different for HIV-infected individuals 
on ART [6) or not on ART (p). 
doi:10.1371/journal.pone.0090985.g001 
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Multiplying the number of susceptible MSM by this transmission 
rate gives the frequency of new infections. Since MSM in Beijing 
often have multiple sexual partners, we classified them as regular 
partners (RP) and casual partners (CP), because they typically have 
different risk profiles [32,35]. We also classified sexual behaviors of 
MSM as insertive or receptive anal intercourse (AI) with difiFerent 
HIV transmission probabilities per sex act [8] . We investigated the 
relationship between per-act and per-partner AI transmission 
probabilities for a given number of sex acts using the Bernoulli 
process assuming independence of risk for each sex act within a 
partnership [38,39]. Parameter notation and descriptions are 
provided in Table 1; initial conditions and parameter values are 
described later. 

Using the Bernoulli process model of HIV transmission during a 
given number of sex acts, the non-infection probability by infected 
regular sexual partners through receptive AI ((1— r)ap and 
insertive AI irdj) under condom-use rate rfj can be described as: 

(1) 

where r is the proportion of sex acts which are insertive AI and 6 is 
the single act transmission probability, and «J and rfj are the mean 
number of sex acts and condom use rates with superscripts r 
(above) or c (not shown) denoting regular or casual sex 
partnerships, respectively. 

The HIV transmission rate is defined as: 

A(o=i-j8;(o-j6;« 

(2) 

where ^{t) is the effective contact rate for all MSM jVat time i with 
Nit) 

The definitions for parameters in the equation (3) are described 
in Table 1. 

From this model, we calculated the basic reproduction number 
(Rq) [40] . Rq is the number of secondary cases produced by a 
typical HIV-infected MSM during his entire period of infectious- 
ness in a demographically steady susceptible population. Calcu- 
lating Rq is critical to determine whether HIV will increase, 
stabilize, or decline among the MSM population in Beijing. 

Data source and estimation 

To simulate the HIV prevalence rates during 2000-2010 
among MSM in Beijing, we used data from local HIV surveillance 
systems and published studies [31,32,33,35,41]. As these data were 
summary statistics, and did not include any identifiers that could 
link the data to individual subjects in the local HIV sur\'eillance 
systems or participants from the published studies, consent was 
waived and the mathematical modeling protocol was approved by 
the institutional review boards of the National Center for AIDS/ 
STD Control and Prevention of Chinese Center for Disease 
Control and Prevention and Vanderbilt University. 

We employed a Metropolis-Hastings algorithm to carry out 
extensive Markov-chain Monte-Carlo simulations for estimating 



the mean values of some unknown parameters [42], including 
yearly recruitment rate (or number of new members into the pool 
of the study population each year), risk-reduction rate of healthy 
MSM (during 2000-2010), relative infectiousness of risk-reduction 
to non-testing subgroups, the per-act probabilit)' of unprotected 
insertive or receptive anal intercourse resulting in an infection, and 
the relative infectiousness ratio of HIV positive MSM in the 
infection stages 2 to 1 (see Figure 1). The algorithm ran for 
1,000,000 iterations, and we adapted the proposal distribution 
after 500,000 iterations using Geweke's method to assess 
com'ergence [42]. 

Initial conditions. The earUest HIV prevalence data among 
MSM in Beijing were available for 2000. The observed HIV 
prevalence rates among MSM in Beijing during 2000-2010 were: 
1.2% in 2000 [41], 3.2% in 2005 [32], 4.8% in 2006 [32], 5.1% in 
2007 [32], 6.5% in 2008 [31,33], 6.8% in 2009 [35], and 7.8% in 
2010 [43]. In order to set the initial conditions, we had to estimate 
the target population of MSM living in Beijing in 2000. According 
the national census in 2000, there were 14 million of people hving 
in Beijing. We assumed half of them were male and 3% of males in 
Beijing in 2000 were MSM, of which half lived in the city. The 
testing rate in 2010 was very low (assumed 0.125%). Also, we 
assumed there was no risk reduction intervention in 2000, as the 
burgeoning HIV epidemic in MSM did not receive attention from 
Chinese public health programs at that time [44]. For these non- 
testing HIV-positive MSM, we assumed one-third were in 
infection stage 1 (CD4 count&350/|J,L) and two-thirds were in 
the infection stage 2 (CD4 count <350/ |J,L). Also, we assumed that 
there was no ART available for MSM in 2000. The resulting 
parameter values for the initial conditions are shown in Table 1 . 

Progression. With an 1 8 year average age of sexual debut, 
we considered a natural removal rate per year which corresponds 
to an a\erage of 42 years of sexual activity and non-AIDS 
mortality (/x = 0.024). The recruitment rate per year into HIV- 
negative MSM population was estimated to range from 2-6.4% of 
the total population [8] giving the recruitment rate (a = 5012) 
(Table 1). The progression duration from HIV infection to AIDS is 
10-12 years in the absence of treatment and may increase by- 
additional 8-12 years with ART [5,45,46]. We used estimated 
disease progression rates being an average duration of 5.6 years 
without ART (therefore, p, = = 0- 1 8), and 11.2 years with ART 
(^1 = (52 =0.09) [5]. We also performed sensitivity analyses for a 
range of disease progression rates. 

Sex and condom use. We estimated condom use rate and 
frequency of AI based on three cohort studies among MSM in 
Beijing [31,32,47] and one meta-analysis among Chinese MSM 
[34]. Table 1 shows condom use rates for insertive and rec('ptive 
AI {rij and rjj) ranging from 30.7%-41.4% [34] and thi; a\-c-rage 
frequency of AI with regular and casual sexual partners (aj and a?) 
ranging from 18.2-36.6 [31,32,47]. 

Transmission. Relative infectiousness of HIV-infected 
MSM (MSM+) who receive HIV testing versus non-testing was 
defined as the ratio of high risk behavior of testing MSM_|_ to that 
of non-testing MSM+. The rationale is that HIV testing (and 
counseling) may lead to change of risky behaviors, such as 
reduction in sexual encounters or increasing condom use, or both. 
Based on risky behaviors among MSM+ from three cohort studies 
[31,32,42], we calculated the relative infectiousness ratio 
(i/'i=0.74; 13.8/18.6) (Table 1). We assumed that non-testing 
MSM did not have different risk behaviors by HIV status (as they 
are unaware of their status) and simulated the effect of risk 
reduction intervention to non-testing MSM (i/'2 = 0.48). The 
conservative estimate for the relative infectiousness among 
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Table 1 . The notation, description, value and uncertainty range of parameters and their data sources for the mathematical model. 







Parameter 


Description 


Value 


Uncertainty 
range 


Source 


Initial conditions 


Sn 


No. of non-tested susceptible MSM 


102,271 








Sf 


No. of tested susceptible MSM 


128 








Sr 


No. of risk reduction (RR) susceptible MSM 


0 








4 


No. of non-tested MSM+ in infection stage 1 


828 








4 


No. of non-tested MSM+ in infection stage 2 


414 








4,4 


No. of tested MSM+ 


1 








4,4 


No. of RR MSM+ 


0 








4.4 


No. of MSM+ on antiretroviral therapy (ART) 


0 






Progression 


a 


Recruitment rate into target population per year 


5,012 


2,073-6,633 


Fit 






Natural removal rate of MSM- per year 


0.024 


_ 


[8] 




PuPi 


Progression rate of MSM+ without ART 


0.18 




[5] 




Si,S2 


Progression rate of MSM+ with ART 


0.089 


- 


[5] 


Sex and condoms 


"n 


No. of anal Intercourse (Al) for non-tested MSM+ with RP 


36.5 


32.2-40.1 


[31,32,47] 






No. of Al for tested/RR HIV_ with RP 


36.2 


34.9-38.6 


[31,32,47] 




<^N 


No. of Al for non-tested HIV_ with CP 


18.2 


11.3-21.7 


[31,32,47] 






No. of Al for tested/RR HIV_ with CP 


22.2 


20.9-24.6 


[31,32,47] 




ffN 


Condom use among non-tested MSM with RP 


30.7% 


28.0-36.7% 


[31,32,34,47] 




I't'Ir 


Condom use among tested/RR MSM with RP 


31.6% 


23.8-37.2% 


[31,32,34,47] 




'fif 


Condom use among non-tested MSM with CP 


37.7% 


34.5-42.3% 


[31,32,34,47] 




lT''fR 


Condom use among tested/RR MSM with CP 


41 .4% 


37.0-43.7%! 


[31,32,34,47] 


Transmission 


'Al 


Relative infectiousness (Rl) of tested MSM+ to non-tested MSM+ 


0.74 


0-1 


[31,32] 




ll'l 


Rl of RR MSM+ to non-tested MSM+ 


0.48 


0-1 


Fit 




e 


Rl of MSM+ on ART to not on ART 


0.4 


0.01 -.4 


[48] 




CO 


Rl of MSM+ in infection stage 2 to stage 1 


1.63 


0.5-3 


Fit 




B,ec 


New infections due to unprotected receptive Al 


0.016 


0.008-0.028 


Fit 




Bins 


New infections due to unprotected Insertive Al 


0.0012 


0.0006-0.0019 


Fit 




r 


Sex acts of insertive Al 


0.56 


0.50-0.61 


Fit 


Intervention uptake 




Testing rates among MSM_ or MSM+ 


0.09 


0.05-0.20 


[49] 




h 


Linl(age to risk reduction (RR) rate among MSM_ 


0.19 


0.10-0.25 


Fit 




l2,h 


Linkage to risk reduction (RR) rates among MSM+ 


0.33 


0.25-0.5 


[49] 




1\ 


ART initiation rate among linked MSM+ in their infection stage 1 


0.018 


0.005-0.036 


[49] 




11 


ART initiation rate among linked MSM+ in their infection stage 2 


0.53 


0.44-0.62 


[49] 




Ml 


ART dropout rate 


0.024 




[49] 




h 


Immunological recovery on ART 


0.50 


0-1 


Fit 



Note: MSM+: HIV-infected MSM; MSM+: 
doi:l 0.1 371/journal.pone.0090985.t001 



HIV-unlnfected MSM; RP: regular sexual partner; CP: casual sexual partner. 



MSM+ who are on ART versus iion-ART (fi = 0.4) was based on 
evidence that ART can reduce HIV transmission by 60 % [48] . 

The per-act risk of HIV transmission through unprotected 
receptive Al ranged from 0.008—0.028, while the risk through 
insertive Al ranged from 0.0006-0.0019 [8]. The estimated 
proportion of practicing receptive anal sex among MSM ranged 
from 0.4—0.5 while the proportion of practicing insertive sex 
ranged from 0.5-0.6 [8], and we estimated the proportion of 
insertive versus receptive sex acts as 0.56 (Table 1). 

Intervention uptake. The average testing rate for MSM 
during 2001-2010 was estimated to be 9% [49]. The linkage rates 
for HIV-positive MSM ranged from 17% to 85% with a mean rate 
of 33% [49]. The ART initiation rates were estimated to be 1.8% 



and 53% among MSM whose CD4 counts are 
<350/|iL, respectively [49]. 



=350/nL and 



Hypothetical scenarios 

We used the model to project HIV prevalence and incidence 
trends during 2011-2020 under five different scenarios: 

S 1. Maintaining the coverage of interventions as they were (i.e., 
average coverage during 2000-2010): ART coverage 
(ti =0.018 and T2 = 0-53), condom use {t]j = rfg^ = Q.'il and 
rf'P = rf^ = QA\), HIV testing (^y = 0.09), and linkage to care 
(h =0.19,/2 = /3 = 0.33) [31,32,34,47,49]. 

S2. Increased ART coverage only {z\ =0.25 and T2 = 0.75). 
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5 3. Increased TLC/ART, but no increased condom use 
(ti = 0.25,T2 = Q.15,ki = 0.25,h = Q.S,h = h= 0.75). 

54. Increased condom use only = ^;5j = 0.63 and 
^5, = ^5^ = 0.83). 

55. Increased TLC/ART (n =0.25,T2 = 0.75,fc, =0.25,/i =0.5, 
72 = 4 =0.75) and condom use (jj^ = »;^ = 0.63 and 
^5, = ^5^ = 0.83). 

In S4— S5, we assumed that condom use doubled since 2000. 

Sensitivity analysis 

Finally, we performed sensitivit}' analyses to assess the 
robustness of the model results. First, by perturbing individual 
parameters, we investigated their influence on estimates of R^. 
Second, we estimated uncertainty ranges for model parameters 
based on literature reviews (Table 1) and randomly sampled 
20,000 parameter sets over these uncertainty ranges using a Latin 
Hypercube Samphng design [50,51] to derive a frequency 
distribution for R(f. Analyses were performed using Matlab 
computer software (MATLAB 7.7, The MathWorks Inc., Natick, 
MA, 2008); analysis scripts and supplemental materials are 
available online (http://biostat.mc.vanderbilt.edu/ 

ArchivedAnalyses) . 

Results 

Using parameter values summarized in Table 1, our model 
produces a best-fit curve to historical HIV prevalence rates among 
MSM in Beijing during 2000-2010 (Figure 2A). The basic 
reproduction number in 2010 is i?o = 2.09 (95% confidence 
interval (CI), 1.83-2.35). Because J?o>l> if the epidemic 
assumptions remain unchanged from 2010, the HIV prevalence 
is expected to substantially increase among MSM in Beijing. Using 
this best-fit model of HIV transmission during 2000-2010, we 
project HIV prevalence for the next ten years (until 2020) under 
five different scenarios. The HIV prevalence is 7.8% in 2010. If 
there is no change in the coverage of interventions or maintaining 
the average level during 2000-2010 (Figure 2, line 1), then the 
HIV prevalence is estimated to increase to 21.4% (95%CI, 9.9- 
31.7%) (Figure 2A, line 1), with an incidence rate of 3.1% (95%CI, 
1.0-5.1%) by 2020 (Figure 2B, line 1). By increasing the coverage 
of ART among MSM^. who are linked to risk reduction (Figure 2, 
line 2), we expect a modest improvement in lowering HIV 
prevalence (19.9%; 95%CI, 9.9-28.4%) and incidence (2.6%; 
95%CI, 0.8-4.6%). If TLC and ART coverage increases (Figure 2, 
line 3), the projected HIV prevalence and incidence will be 14.5'X> 
(95%CI, 7.0-23.8%) and 1.4% (95%CI, 0.2-3.3%). If MSM 
increase condom use only (Figure 2, line 4), then HIV prevalence 
and incidence will be 13.0% (95%CI, 9.8-28.4%) and 1.5% 
(95%CI, 0.8-4.6%), respectively. Bundling increased TLC/ART 
with increased condom use (Figure 2, line 5) may result in an HIV 
prevalence of 8.7%. (95%.CI, 5.4-1 1.5%) and an incidence of 0.5% 
(95%CI, 0.1-0.9%) by 2020. It is important to note diat only this 
bundled HIV intervention strategy is likely to have an 
J?0 = 0.93<1 and eventually result in HIV eradication in this 
population. 

We performed simple sensiti\'ity analysc's for intervention and 
transmission parameters. First, we analyzed the influence of the six 
intervention parameters (testing and linkage to risk reduction) on 
the prevalence during 2000-2020 (Table SI). The parameters 
related to HIV testing (fc) have the greatest influence on the 
estimates of i^o- In particular, the partial rank correlation 
coelficient (PRCC) values range from —0.98 to —0.59, indicating 
that 7?o decreases for higher testing rates. Similarly, the parameters 



related to HIV linkage to risk reduction (i) also impact Rq. PRCC 
ranged from —0.44 to —0.14, indicating that Rq decreases for a 
higher hnkage rate. Figure 3 shows two contour plots of J?o under 
diflerent HIV testing or linkage to risk reduction among MSM+ 
while condom use remains at the average level during 2000—2010. 
Even if all MSM+ are tested (testing rate = 100%) or are linked to 
risk reduction (RR rate = 100%) without increasing condom use, 
the reproduction number remains above the epidemic threshold 
(^o>l). 

Second, we analyzed the influence of the relative infectiousness 
of MSM+ on ART versus non-ART (varying 0.01-0.4) and that of 
MSM+ with CD4 counts <350/nL versus >350/|iL (varying 0.5- 
5.0). In all perturbations for relative infectiousness, Rq> \ endures. 

Third, we found that Ro>\ is sustained if MSM increase 
condom use but the coverage of TLC interventions remains at the 
a\-erage level during 2000-2010. 

Fourth, we varied the average duration of disease progression 
without ART from 3.7-1 1.1 years and with ART from 7.4-22.2 
years. The median J?o is 2.1 1 and consistently greater than 1. 

In extensive sensitivity analyses, 20,000 parameter sets were 
randomly sampled from corresponding uncertainty ranges (Table 
S2). The distribution of R^ is skewed to the right with mean 1 .93 
(standard deviation: 0.82) and median 1.82 (interquartile range: 
1.36-2.42). A small percentage (1 1.4%) of estimates falls below the 
epidemic threshold (J^o < 1). Relative infectiousness of RR to non- 
testing MSM+, transmission probability of unprotected receptive 
AI, and relative infectiousness of MSM+ on ART to non-ART are 
the three most influential parsimeters with PRCC equaling to 0.95, 
0.90 and 0.71, respectively. 

Discussion 

Combination testing and treatment inter\'ention packages have 
been considered as a promising strategy' for HIV prevention 
among MSM in recent years [52]. We explored the effects of 
TLC/ART, along with condom use, for the emerging HIV 
epidemic among Chinese MSM [4,53]. Our projections predicted 
that with the current level of services HIV prevalence among 
MSM in Beijing will almost triple from 7.8% in 2010 to 21.4% by 
2020. Expansion of ART coverage in the current continuum of 
care among MSM+ would only lead to a moderate reduction of 
this rising trend, though enhanced TLC interventions plus 
increased ART coverage will further reduce incident HIV 
infections. Only with a combination of increased condom use 
and a higher TLC/ART coverage might HIV epidemic be 
expected to decline. Hence, promoting safer sex education and risk 
reduction counseling service to increase condom use is crucial for 
HIV prevention among Chinese MSM, even in the context of 
expanded TLC and ART. 

As obser\'ational studies and the HPTN 052 randomized 
controlled trial have demonstrated the effectiveness of ART for 
reducing heterosexual HIV transmission [54,55,56,57,58], treat- 
ment as prevention has been advocated as a strategy for HIV 
prevention among MSM [59] . However, the impact may vary by 
epidemic scenarios [21]. In settings with a long-standing epidemic 
like San Francisco, TLC with ART expansion may reduce the 
HIV epidemic [52]. Our model is less optimistic: that ART 
interventions alone will succeed in substantial prevention reduc- 
tion. The HIV epidemic among Chinese MSM is relatively recent 
[44], and the majority of cases may stiU be in earlier stages of 
disease, and therefore have not received ART [13]. We use local 
public health data such as ART coverage rates for estimating the 
parameters in the mathematical models; therefore our modeling 
results likely reflect what is actually happening in Beijing. A 
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comparable example is that in the HPTN 052 trial earlier ART 
treatment reduced heterosexual transmission by 96% [7], while 
among 38,862 Chinese serodiscordant couples the rate of 
transmission was only reduced by 26% in the first year of 
treatment and the reduction lost statistical significance in 
subsequent years [60] . The former result comes from a controlled 
trial setting while the later one is from a real-world public health 
setting. Our model projections are consistent with the study in the 
United Kingdom where increased condom use was more likely to 
drive down HIV incidence than high ART coverage [20] . 

Increasing coverage of HIV testing may lead MSM+ in Beijing 
to reduce risky behaviors, notably the use of condoms; partner 
reduction interventions and serosorting (i.e., MSM of the same 
HIV status) interventions could conceivably have a similar benefit, 
but we did not model these factors. Based on our model, HIV 
prevalence will continue to increase, as fewer MSM+ die and fall 
out of transmission networks. The impact of HIV testing is largely 



dependent on how those infected MSM with known status reduce 
their risk behaviors, particularly practicing anal sex without 
condom use. Given that only a fraction of HIV-infected men know 
their status, are eligible for ART, and successfully negotiate the 
system to adhere to ART, increased condom use is a crucial 
element in the combined interventions for reducing the epidemic. 

The limitation of our mathematical model is a simplification of 
the real world epidemic. Like other models, we balance parsimony 
and reality. For example, we assumed that all linked MSM+ 
receive risk reduction. This is consistent with Chinese HIV testing 
protocol that all MSM+ are expected to receive risk reduction 
counseling. However, in reality, the quality of counseling may vary 
by chnic and client. Because some counseling may not lead to 
behavioral changes, we considered condom use as an independent 
factor in our model. Although model parameters are carefully 
specified using the available literature, they may not reflect the 
current realities among Beijing MSM. For example for scenario 3, 
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HIV-infected Chinese MSM while condom use remains at the average level of years 2000-2010. (A) Under the different coverage rates of 
HIV testing intervention; (B) Under the different coverage rates of HIV risk reduction (RR) intervention. 
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we assumed increasing HIV testing from 9% to 25%; increasing 
the linkage rate to risk reduction education from 19% to 50% 
among HIV-uninfected MSM and from 33% to 75% among HIV- 
infected MSM; increasing ART coverage from 1.8% to 25% 
among HIV-infected MSM whose CD4+ cell count a350/|J,L and 
from 53% to 75% among those with CD4 count <350/|iL. The 
magnitude of these assumed increases in coverage are based on 
our best estimations of how well the public health system could do, 
in the view of local HIV/ AIDS experts (including authors YR, 
XH, HL and YS as well as non-author experts from local centers 
for disease control). However, we performed extensive sensitivity 
analyses by varying input parameters and found that our findings 
are robust to reasonable perturbations. 

No scenario can result in a substantial decline in HIV epidemic 
by 2020 among Chinese MSM, except the bundled package of 
transmission prevention strategies including increased TLC/ART 
and condom use. Simulations show that TLC alone has some 
impact for HIV reduction, but far less than increasing condom 



use. Our model suggests that increasing condom use must be a 
component of combination intervention packages to achieve 
significant reduction in HIV incidence. 
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